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We analyse the linear spin wave theory calculation of the superfluid phase of a hard-core boson 
model with nearest neighbour exchange interaction J, Dzyaloshinskii-Moriya (DM) exchange inter- 
action D and four-particle ring-exchange interaction K at half filling on the triangular lattice, as 
well as the phase diagrams of the system at zero and finite temperatures. We find that the DM 
interaction can be removed from the Hamiltonian by redefining the spin operators but this leads to 
a change in the nearest neighbour exchange interaction. We also find that the pure J model (XY 
model) which has a well known uniform superfluid phase with an ordered parameter (Sf) ^ at 
zero temperature is quickly destroyed by the inclusion of a negative- if ring-exchange interactions 
for D 0, favouring a state with a (^, 0) ordering wavevector. We further study the behaviour of 
the zero temperature superfluid density and finite-temperature Kosterlitz-Thouless phase transition 
(Trt) in the uniform superfluid phase for some values of k = K/ J, 77 = D/ J, by forcing the universal 
quantum jump condition on the finite-temperature spin wave superfluid density. At zero tempera- 
ture, we find that the maximum values of the superfluid density as a function of k increases as r\ 
increases which shows that the DM exchange interaction constant increases the zero temperature 
superfluid density. 



I. INTRODUCTION 

The effective studies of continuum field theories have 
resulted in detailed predictions for the low-energy physics 
of quantum spin systems in two dimensions. Spin wave 
theory can provide us with a rather accurate picture of 
the low-lying states of quantum spin systems. There 
are several versions of spin wave theory. The stan- 
dard spin wave theory is based on Holstein-Primakoff 
representation 10 which was first applied to the study of 
Heisenberg model by Anderson^ and further extended to 
second order by Kubc0 and OguchP. 

Spin wave theory was thought to be unsatisfactory 
in the case of XY model until Gomez-Santos and 
Joannopoulous^ showed that by a good choice of quan- 
tization axis, one can obtain a good theoretical result 
for XY model. Since then numerous applications of spin 
wave theory have been carried out on XY model with 
different lattice configurations and the results obtained 
so far are in a good agreement with quantum Monte 
Carlo simulations (QMC)P^. Magnetic spin models can 
be supplemented with antisymmetric spin exchange in- 
teraction called Dzyaloshinskii-Moriya (DM) interaction 
which arises from spin-orbit coupling. This interaction 
was first proposed by DzyaloshinskiP^ and later derived 
by MoriysP^. We shall see later that such additional in- 
teraction does not change the form of the original Hamil- 
tonian, it only leads to changes in the exchange interac- 
tion constants. 

Another interesting area is the multiple (ring) spin- 
exchange models. This model was first introduced to 
describe the magnetic properties of solid 3 He^. It incor- 
porates ring-exchange interactions over plaquettes such 
as spin 1/2 four-spin XY ring exchange of the form: 



H K = -K^2 (S+SrS+Sf + SrS+S^S+) , (1) 

(ijkl) 

where the summation runs over plaquettes with the in- 
dices running counter-clockwise, i and j are nearest 
neighbours lying opposite to k and /. 

The ring exchange interaction is important in Wigner 
crystal near the melting density 13 . This model, alone or 
in competition with pure XY model with nearest neigh- 
bour exchange has attracted considerable attention over 
the years and has been studied extensively on a square 
lattice using a stochastic series expansion (SSE) quan- 
tum Monte Carlo method 6 also a comprehensive theoret- 
ical study (spin wave theory) has been done on a square 
lattice 8 . It has been suggested in recent works that mod- 
els of this form may harbour exotic ground state prop- 
erties, including de-confined quantum critical points or 
quantum spin liquid phases. 

In this paper, we shall calculate the exact (numeri- 
cal) value of the Kosterlitz-Thouless temperature, super- 
fluid density and other low temperature thermodynami- 
cal properties of spin 1/2 XY model with ring exchange 
and DM interactions on a triangular lattice using linear 
spin wave theory without the inclusion of the repulsive 
interaction between bosons. 

The format of the paper is as follows: In Sec. II, we 
present the model Hamiltonian. In Sec. Ill, we apply lin- 
ear wave theory by choosing our quantization axis along 
the ^-direction and use it to diagonalize the Hamiltonian 
and obtain its energy. In Sec. IV, we analyze the disper- 
sion and plot it for some values of k and 77. In Sec.V, 
we explore the zero temperature superfluid density using 
the diagonalized Hamiltonian and plot it as a function of 
k for some values of 77 . In Sec. VI, we calculate the finite 
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temperature superfluid density for some values of k and 
r] and the values of the Kosterlitz-Thouless temperature 
for this model and in Sec. VII, we make some concluding 
remarks. 



II. MODEL 

In this section, we shall present the model Hamiltonian 
and the mean field theory argument of the J-K model on 
a triangular lattice, similar work was done on a square 
lattice 8 for D = 0. Our model Hamiltonian is given by 

(ij) 

(ijkl) 



H 



-K 



-D-iSixSjf) 



(2) 



where & = (Sf, S?) and Sf = S? ± zSj. 

The first summation is over nearest neighbour pairs on 
a triangular lattice and the second summation runs over 
the three possible plaquette orientations on a triangular 
lattice as shown in Fig. 0. It has been shown that this 
Hamiltonian undergoes a Kosterlitz-Thouless phase tran- 
sition for Z), K = at Tkt ~ 0.69 for 2D model and a 
superfluid phase for temperatures less than Ttk 16 - For 
D = 0, K < 0, there is a sign problem which prevents 
(QMC) simulations^ though it is not possible to capture 
this sign problem in the linear spin wave theory. QMC 
simulation on a square lattice has been performed for 
D = 0, K > in which there is no sign problem^- For 
simplicity, let's assume that D = De z , that is D points 
in the ^-direction. Then ([2| takes the form 



H 



^[(J + iD)StSr- 

(ij) 

-K (StS-S^Sf 

{ijkl) 



(J-iD)S~S+] 



(3) 



Rotating globally about the z-axis^by an angle ±a/2 
we have 



/2 o+ 



S + e -ia/2 



si 



S+e~ ia / 2 ,S+ -^S+e ia/2 . 



(4) 



where a — arctan(Z)/ J). Plugging Q into ([3|, the 
Hamiltonian transforms into 



-K ^2 (S^Sj ^k^i + $i Sj~$k ^ + ) • 

(ijkl) 



(5) 



where J D = V J 2 + D 2 = Jyfl + rf, rj = D/J. 

Therefore we see that the redefinition of the spin oper- 
ators removes the DM interaction and rescales the near- 
est neighbour interaction but has no effect on the ring 
exchange term. 



Lets define the spins as classical vectors by making 
the transformation = pe 1 ^*. The model Hamiltonian 
becomes 



H = -2J D ^2 P 2 cos ~ fa) 

(ij) 

- 2K ^ P 4 cos (fa _ fa + fa ~ fa) • 

(ijkl) 



(6) 



Now lets consider the case Jd-, K > 0, in this case, min- 
imizing the energy we have fa = <fij for the Jd term and 
fa = 4>j, fa = fa for the K term which leads to a ferro- 
magnetic ordering of spins. Consider the case J^,i^<0, 
in this case, minimizing the energy we have fa — (j)j = tt or 
fa = 0, fa = 7r for the Jd term and fa — + fa — fa = 
or fa = fa = fa = 0, fa = 7r for the K term which leads 
to ud (up-down) state for the Jd term and uuud state 
on the plaquettes, basically we have two configurations 
of spins on the lattice. 

This model has also been considered by L. Balents and 
A. ParamekantP^ with the inclusion of the repulsive in- 
teraction U term between bosons without the DM term 
D = 0, that is Jd = J and also in the regime where 
J « K. They showed that when J = 0, the four-spin 
exchange leads to a manifold of ground states with gap- 
less excitations and critical power-law correlations and 
when J ^ 0, fluctuations select a four- fold ferrimagneti- 
cally ordered ground state with a small spin (superfluid) 
stiffness which breaks the global U(l) and translational 
symmetry but they did not obtain any exact (numerical) 
value of the Kosterlitz-Thouless temperature and the su- 
perfluid density. 

From their Hamiltonian they argued that with J = 0, 
the ground state is independent of the sign of K and for 
non-zero J, the sign of K is vital, for J < with K < 0, 
there is a ferromagnetic phase while J > leads to a 
y/3 x y/3 Neel order which are also the ground states for 
large \J/K\. They finally concluded that there are no 
phase transitions at any non-zero J other than the well 
known phases. 



III. LINEAR SPIN WAVE THEORY 

The basic assumption of spin wave theory lies on se- 
lecting a classical ground state and determining the fluc- 
tuation around it. In other words, one considers quantum 
fluctuations very close to an ordered ground state config- 
uration of the system under study. By the usual mapping 
between spins 5=1/2 and the hard-core bosons, we can 
view Q as a hard-core boson model. For J » D,K or 
D,K = 0, the T = ground state (in-plane ferromag- 
net) has an ordered parameter M x = (Sf) ^ which 
breaks the U(l) global rotational symmetry or a super- 
fluid phase in the hard-core boson versiorP^^. One can 
therefore perform a spin wave expansion around this or- 
dered state configuration by introducing the boson oper- 
ators CLi and a\ which represent the low-energy spin wave 
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FIG. 1: (Color online) The three plaquette orientations (a), 
(6), and (c) on a triangular lattice with the position coordi- 
nates S a = (1,0), tip = ^ 7 1 
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excitations out of (Sf) and treat other terms in ([5| as 
perturbations. 

We shall follow the procedure of Gomez-Santos and 
Joannopoulos^ and choose our quantization axis along 
the x direction (instead of z direction). This allows us 
to write the Holstein-Primakoff representation 10 in the 
linear spin wave theory as 



Sf — - cl[clu 



(7) 



The method for diagonalizing the Hamiltonian (5l) is 
given by R. Schaffer et al 8 . It involves writing (5l in 
terms of Sj and Sj using Sf = ± iS?, after which 
we substitute ^ into the resulting equation and then 
Fourier transform taking into account the summation 
over all the three plaquettes for the if -term and finally 
diagonalize using Bogoliubov transformation. After all 
these mathematical gymnastics, the diagonalized Hamil- 
tonian is 



H = H MF + Yl ~ Ak ) 

k 



(8) 



Here, the mean-field energy and the coefficients are 
totally different from those obtained on a square lattice^, 



they are given by 
Hmf = 



3[ ±J D N+^KN 



Al-B* 



where 



A k = J D Qk + KR k , 
^k = JdSil + KT^. 



Q k = 3 (l - ^ 

Q 3 

1 1 1/ 

2 - 2^k + g (7k + 7k) J- , 



r k = 3<ji 7 k-^(7k + 7k) 



(9) 

(10) 

(11) 
(12) 



(13) 
(14) 

(15) 

(16) 



and the lattice structure constants are given by 



?k = - I cos k x 



7k 



- | cos V3k y 



k x y/3k y \ 
2 cos — cos , 

2 2 y 

3k x ^/3k y 

2 cos cos 

2 2 



(17) 



We can see that the Hamiltonian (|8| reduces to pure XY 
model in the limit K = as expected. This diag- 
onalized Hamiltonian ([8| will be used to analyse some 
properties such as dispersion, ground state, internal en- 
ergy, superfluid density etc., all as a function of n and 

7]. 



IV. SPIN WAVE ENERGY DISPERSION 

We shall now proceed by analysing the dispersion of 
this model for some values of n and 77. We shall set J = 
1/2 which corresponds to the parameter value of pure 
XY model for which D, K = 0. From the diagonalised 
form of the Hamiltonian, the dispersion for n = —4/3 is 
given by 



c (k) = 2cj k = 2 



7k- 



(18) 



The graph of ( p~8| ) is shown in Fig.Q for 77 = 0. It shows 
a zero mode at the corners and the center of the Brillouin 
zone of a triangular lattice. The dispersion along the k x 
direction for k y = 0, 2tt/^ and k = —4/3 is also shown 
in Fig.Q for several values of 77. There are three zero 
modes at k=Q=(0,0), k= ±Q = (±^,0) for 77 = 0, k y = 
0, away from this value of 77, there is only one zero mode 
which is at the center of the Brillouin zone k=Q=(0,0). 
For 77 = 0, k y = 27r/v // 3, there is only one zero mode at k= 
±Q =(±27r/3,0). If one chooses ferromagnetic ordering 
along the k y axis then the corresponding ordering wave 
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vector is k= ±Q = (±4^,0) which is the soft modes of 
the dispersion for ft = — 4/3, r] = 0. The linear spin wave 
instability of the excitation spectrum at the corner of the 
Brillioun zone k= ±Q = (±^,0) occurs for 77 = 0, k = 
—4/3. For the pure XY model 77 = 0, K/J = 0, there is 
a gapless excitation at Q =(0, 0) (Goldstone mode of the 
superfluid phase) , but there is no minima at the ordering 
wave vector. Near the zero modes for 77 = 0, n = —4/3 
we have 



e(k) 



{kl + k 



2\V2 



V3Jk. 



(19) 






FIG. 2: (Color online) The dispersion e (k) for 77 = and k — 
—4/3 (first). The dispersion along k x direction for several 
values of 77 for k = —4/3, k y = 27r/v / 3 (second), and for 
k = -4/3, k y = (third). 



V. SUPERFLUID DENSITY AT T = 

In this section we shall explore the superfluid density 
of this model on a triangular lattice. This is done by 
making a phase twist on the lattice sites. The rele- 
vant transformations have been worked out for a similar 
model on the square lattice^, with this transformations it 
was found that the ring exchange term is independent of 
the phase angle. Using this transformation, the twisted 
Hamiltonian becomes 



H{9) = -2J D { (SfSj + S^V) cos 

(ij) 



■(SfS? -S?S%) sm9} 
K (S+S-S+S-+S-S+S^S+) 



(20) 



(ijkl) 



where = 0i — 0j. In the linear spin wave theory, the 
coefficient of sin does not contribute, therefore we have 



H(0) = -2J D (SfS* + S?S?) cos (9 



(ij) 



(21) 



(ijki) 



Comparing ([5| and (21 ), we see that the effect of the twist 
is to rescale the nearest-neighbour exchange interaction 
in ([2]) as Jd — » Jd cos 0. Therefore, the diagonalized 
form of (21) is exactly of the form (pi) with Jjj replaced 



by Jd cos 0. That is 

H{6) = H MF (0) 



k 



K(0) - A k (6)) 



^2 Wk ( 6 ') ( a lc a k + al k a-k) 



(22) 



where the twisted mean-field energy and the coefficients 
are given by 



Hmf(8) 



-3 ( ^J D Ncos8 



1 



-KN 



A k (6) = J D Qk cos 6 + KR k , 
B k (6) = J D S k cos6 + KT k , 



(23) 



(24) 



(25) 



the coefficients Q k , R k , Sk, and T k remain the same as 
dll-p). 



(26) 



The free energy is given by 
F(0) = ~\nZ(0)=Eo{e) + ±hi (l 
where the ground state energy E (6) is given by 

E O (0) = H MF (e) + k(^) - M^)) . (27) 
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At zero temperature, the free energy is simply the ground 
state energy. Taylor expanding the ground state energy 
we have 



E (9) = E (6 = 0) + \p s 2 + O{0 A ), 



Ps(T = 0) 



1 d 2 E (6 ) , 

n oe 2 



6>=0 • 



(28) 



(29) 



The expression for p s can be obtained from ( |27| ) by 
expanding the cos to order 2 and differentiating. This 
gives 



Ps(T = 0) = 



Jd 
2N 



i 



(30) 



(Jd(Q 



2 si) 



K(Q k R k - S k T k )) 



where we have divided by 2 to account for the dimen- 
sionality of the lattice. 

The plot of p s (T = 0) is shown in Fig.(3| for a range 
of k and several values of rj. The superfluid density curve 
shows a maximum at k = 0, with a value of p s = 0.4059 
for 77 = 0, p s = 0.4140 for 77 = 0.2, p s = 0.4372 for 
77 = 0.4, p s = 0.4734 for 77 = 0.6. We see that the max- 
imum of the superfluid density increases as r] increases, 
therefore the effect of the DM interaction is to increase 
the superfluid density. We also notice that the curves 
decrease monotonically with increasing \K\ as one moves 
away from these maxima, in other words the ring ex- 
change term decreases the value of p s (T = 0) from the 
pure XY result. On the positive- K side, p s decreases rel- 
atively gradually, only becomes zero for extremely large 
values. This is consistent with the result obtained in the 
dispersion, which indicates that no soft modes develop 
for moderate values of positive K. Similar result was ob- 
served on a square lattice 8 for r] = 0, which is consistent 
with quantum Monte Carlo simulation. On the negative- 
K side, the value of p s decreases rapidly as it approches 
k = -4/3. 




FIG. 3: (Color online) The superfluid density as a function 
of k for several values of r/ (top). The mean- field value of 
superfluid density is pMF = 0.375 for 77 = 0. 



VI. SUPERFLUID DENSITY AND 
KOSTERLITZ-THOULESS TRANSITION AT 
FINITE TEMPERATURE 

In this section we shall discuss the uniform superfluid 
phase of our model at finite temperatures. In order to 
calculate the finite temperature superfluid density, we 
replace the twisted ground state energy in (29) with the 
twisted free energy, that is 



N dO 2 



Using (26) we obtain 



3Jjj Jd 
~T~ 27V 



(31) 



(32) 



(A k Q k - £ k # k ) 1 



WT _ 1 



This expression for p 3 (T) is plotted as a function of T/J 
in Fig. (|4|) for some values of k and 77. 





FIG. 4: (Color online) The superfluid density as a function 
of T/J for k — (top), lines are labelled by the parameter 
value 77 and for k ^ (bottom), lines are labelled by the 
parameter value 77 and k . The dashed line is the universal 
jump condition. 

The graph shown in Fig.Q is similar to the one ob- 
tained on a square lattice^ for 77 = which shows slowly 
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decaying superfluid density. We see that the superfluid 
phase for k = 0,7? = is being destroyed for k < re- 
gardless of the value of rj. The dash line is the so-called 
universal jump condition 

fi^ST) _ ? (33) 

which accounts for the discontinuity in p s {T). The esti- 
mate of Tkt can be found by solving p s (T)/T = 2/-7T, 
using p s (T) from our spin wave theory. Using this 
procedure we found the KT transition temperature at 
T KT = 0.9295 J for n = 0, rj = and J = 1/2, the pa- 
rameter values for the pure XY model. This is shown 
in Fig.Q, where the dashed line crosses the curve for 
k = 0, r] = 0. We can as well read off the values of Tkt 
for ft ^ 0, r\ 7^ 0. The full analysis of this model for r\ = 
has been worked out previously by the same author of 
this article^. 



VII. CONCLUSION 

We have analysed the study of linear spin wave theory 
of hard-core bosons (zero field XY model) with DM and 
ring exchange interactions at half filling on the triangular 
lattice. We studied through linear spin wave theory the 
destruction of uniform superfluid phase by the bosonic 
ring exchange model on a triangular lattice at half fill- 
ing. The dispersion of this model was calculated by ap- 
plying the traditional Holstein-Primakoff representation 
and summing over the three plaquettes orientations on a 



triangular lattice. This calculation showed a spin wave 
instability and a development of three minima at k= ±Q 
= (±^,0), andk = Q= (0,0) for k = -4/3, rj = 0. One 
should expect a phase transition from a superfluid phase 
to another phase at this wave vector k= Q =(^,0) for 
ft = —4/3,7? = 0. We hope that quantum Monte Carlo 
data should provide further insight into this issue. 

The mean field superfluid density and ground state 
(zero temperature) spin wave superfluid density obtained 
in this model for ft = —4/3, r] = is bigger than that of a 
square lattice 8 despite the fact that both are 2D systems. 
This might be due to a larger number of nearest neigh- 
bours and plaquettes on a triangular lattice. We calcu- 
lated the finite temperature uniform superfluid density 
and use it to estimate the Kosterlitz-Thouless transition 
temperature by forcing it to obey the universal quan- 
tum jump condition. We also analyse the effect of DM 
interaction on the superfluid density by plotting it for 
several values of 77 as a function of ft. It was found that 
the DM interaction increases the maximum of superfluid 
density at zero temperature. It has been shown with the 
model Hamiltonian (|5| on a square lattice using QMC 
simulations for 77 = that for small K>0 away from the 
XY point, the zero-temperature spin stiffness value of 
the XY model is decreasecP. Our results above seem to 
agree with this trend found in QMC simulations. 
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